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Abstract 

O 
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CNJ , The Un — * v T oscillation hypothesis will be tested through v T production of r 

(3J[)! in underground neutrino telescopes as well as long-baseline experiments. We 

provide the full QCD framework for the evaluation of tau neutrino deep inelas- 
tic charged current (CC) cross sections, including next-leading- order (NLO) 
corrections, charm production, tau threshold, and target mass effects in the 
collinear approximation. We investigate the violation of the Albright-Jarlskog 
relations for the structure functions -F4.5 which occur only in heavy lepton (r) 
scattering. Integrated CC cross sections are evaluated naively over the full 
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phase space and with the inclusion of DIS kinematic cuts. Uncertainties in our 
evaluation based on scale dependence, PDF errors and the interplay between 
<^> kinematic and dynamical power corrections are discussed and/or quantified. 
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J3 • I. INTRODUCTION 

a: 

Ph. Results from the SuperKamiokande underground experiment measuring the atmospheric 

neutrino flux suggest that muon neutrinos oscillate into tau neutrinos with nearly maximal 
mixing Jl| . A test of the oscillation hypothesis is v T production of r through charged current 
interactions, a process which will be studied in underground neutrino telescopes as well 
as long-baseline experiments || measuring neutrino fluxes from accelerator sources. For 
precision measurements of oscillation mixing angles and eventually CP violation, neutrino 
cross sections will ideally be know to the level of a few percent. Eventually, measurements 
of neutrino-nucleon charged current interaction cross sections are expected to be at the 1% 
level at a neutrino factory [|]. 

The QCD theory of deep inelastic cross sections in the leading-twist approximation has 
proceeded to the point of evaluating next-to-next-to leading order (2 loop) perturbative 
QCD corrections to the coefficient functions for the structure functions F\, F 2 and _F 3 



X 



and approximations for the splitting functions fl-pl- In the specific case of z/^-isoscalar 
nucleon (JV) cross sections, target mass corrections and nuclear binding effects in a LO 
twist-2 approach of DIS, combined with the elastic peak and a modeling of higher twists 
in the continuum and resonance region of DIS, have also been investigated [12, [L3]. Tau 
neutrino charged current interactions with nucleons have received less theoretical attention 
12, I4],|l5|. Albright and Jarlskog, in Ref. ||16|| , pointed out that there are two additional 



structure functions, F 4 and F 5 that contribute to the tau neutrino cross section. F 4 and F 5 
are ignored in muon neutrino interactions because of a suppression factor depending on the 
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square of the charged lepton mass (me) divided by the nucleoli mass times neutrino energy, 
mf/ (M N E y ). At leading order, in the limit of massless quarks and target hadrons, F 4 and 
F 5 are 

F 4 = (1) 
2xF 5 = F 2 , (2) 

where x is the Bjorken-x variable. These generalizations of the Callan-Gross relation F 2 = 
2xF\ are called the Albright- Jarlskog relations. As with the Callan-Gross relations, the 
Albright- Jarlskog relations are violated from kinematic mass corrections and at NLOQ in 



QCD. We quantify the violation in Section III A 



Below E v ~ 10 GeV for muon neutrinos and higher energies for tau neutrinos, untangling 
the exclusive and inclusive contributions to the neutrino- nucleon cross section is difficult. 
The cross section in this energy range has quasi-elastic, resonant production (such as A 
production), non-resonant pion production and other inelastic contributions. There are 
several phenomenological methods for avoiding double counting in this region. One method 
employs a cutoff in the hadronic invariant mass W > W m i n such that the combined exclusive 
and W m in dependent inclusive cross section yields a total consistent with data |T7[. Using 



muon neutrino data, a common choice for W m i n is 1.4 GeV. A second method involves 
looking at different final state multiplicities fll4fl . By normalizing the calculated total cross 
section for a given multiplicity j, of ot , to the data, one can determine the factors fj in 
°tot = ^rcs + fj a ms- The quantity a 3 ms is determined by the inelastic cross section and the 
hadronization scheme. The Monte Carlo for the Soudan experiment does the normalization 



at E v = 20 GeV [14 



We present here the deep-inelastic contribution to u T N — > rX incorporating next-to- 
leading order QCD corrections, collinear target mass and charmed quark mass corrections. 
The NLO structure functions F 4 and F 5 including charm quark production have been eval- 



uated by Gottschalk in Ref. [|1^]. Here, we correct a misprint and evaluate gluonic helicity 
states in D dimensions rather than 4 dimensions. We show numerical results for the new 
structure functions and for the charged current cross sections, both with and without im- 
posing a nonzero value for W min . To gauge the effects of perturbative delicacies such as 
higher-twist beyond the inclusion of target mass effects in the scaling variable, we also eval- 
uate the effect of imposing a cutoff on Q 2 . We find that the collinear cutoff effect of the tau 
mass m T = 1.78 GeV is enough to guarantee that for E Vt > 5 GeV, <jqq(v t N) is dominated 
by Q 2 > 1 GeV 2 . While this is widely considered a perturbative scale in QCD, some caveats 
about the possible importance of higher twist are discussed. 

In Section O, we show the formulae for the differential cross section in terms of struc- 
ture functions, and the expressions for the structure functions at NLO. Charm quark mass 
corrections for the charged current process are displayed, together with the results for the 



m c — > limit. In Section |TJ, we exhibit our numerical results for the structure functions F4 
and F 5 , and for v T N and v T N interactions for neutrino energies up to 100 GeV. Cross sec- 
tions are compared to those with incident muon neutrinos and uncertainties are discussed. 



1 We will find below that Eq. (0) is not violated in massless NLO QCD. 
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We conclude in Section |V[ In an Appendix we rewrite the Albright- Jarlskog relations in 
terms of helicity amplitudes and pinpoint the approximations in their derivation. 



II. FORMULARY: v T DIS AT NLO 



Neglecting neither the target nucleon mass nor the final state lepton mass m T , the 
charged current u T (anti-)neutrino differential cross section is represented by a standard set 
of 5 structure functions fL6| pl 



dVg G 2 F M N E V 
dxdy ~ tt(1 + WTWf 
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(3) 



where {x, y, Q 2 } are the standard DIS kinematic variables related through Q 2 = 2M N E v xy 
and where we have neglected factors of m 2 /2M N E u ■ Q 2 /M^. These latter corrections from 
the ~ q fJ- q u /M$ v part of the massive boson propagator are negligible both at low and at 
high neutrino energies and will not enter our numerics. For completeness of this formulary, 
though, they can be included multiplicatively by replacing: 



F 



w ± 



F^x 



(4) 



with 



ei 

^3 

e 4 

^5 



w 2 (Q 2 + 2M 2 W ) 
2M^ 
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(5) 



Kinematics determine the integration ranges^] to be [[12], [UJ 



2 Our normalization of F± differs from that in 16] by a factor of x. 

3 A typographic error in Refs. [0,ll6|] is corrected here in the lower limit for x. 
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2M N (E U - m 7 



< x < 1 



(6) 



and 



a — b < y < a + b 



(7) 



where 



1 — m 



\ 



t I 2M N E„x ~ 2El 
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(8) 
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2M N E v x 



2 1 + 



In the perturbative regime we can calculate the structure functions i^ w± from parton 



dynamics by use of the factorization theorem 



(9) 



relating the hadronic and partonic {Cj^ u ) forward matrix element of the product of 

weak currents < J^J U >. In Eq. (|9|), /(£,/i 2 ) is a parton distribution function evaluated at 
factorization scale /i and the parton momentum fraction £ of the light cone momentum of 



the nucleon P£ = (P^ + Pfj)Jy2. Intrinsic transverse momentum of the incoming parton 



is neglected throughout our discussion, i.e. p± — in (J|). 

For the neutrino energies of interest here, we can safely restrict the consideration to the 
first two quark generations. The light flavour contributions to Eq. (||]) can be obtained from 
the m c — > limit of the charm production component which we will, therefore, consider first. 
The charm production contribution to da will be represented by the structure functions F?. 
We introduce theoretical structure functions^ 



^(x,Q 2 ) = (1 - fo) - a^, a* 2 ) + 




(10) 



4 To meet with an experimentalist's point of view, we denote functional dependence on x and Q 
only in the LHS of Eq. (jlOj) ; corresponding in theory to an exact knowledge of all mass parameters 
and an ideal validity of the renormalization group for observable cross sections. Residual scale (//) 
dependence will be investigated below. 
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for scattering off the CKM-rotated weak eigenstate 

s' = \V s f s+\V d f d (ff) 

and its QCD evolution partner 

g'=g(\V s>c \ 2 + \V dtC \ 2 ) (12) 



i.e., ds'/d\nQ 2 = s' ® P qq + g' ® P, 



19 ■ 



In Eq. ( |i0|) we set the renormalization scale equal to the factorization scale and 



77 II 1 M 2 

^ = A '^ = 2^ + V^ + ^ (13) 



is the target mass corrected slow rescaling variable. The quantity rj is the Nachtmann variable 
fi~9fl . The charm mass dependence is included in the dimensionless A = Q 2 /{Q 2 + rn 2 ). In 
Eq. ( |T0D the convolution variable £ in Eq. (^j has been traded for £' = fj/£ which relates to 
the partonic CMS energy s = (p + q) 2 through l/£' = A(l + s/Q 2 ). 

The theoretical structure functions in Eq. (|10D are obtained from a tensor projection 
on the partonic Cj^ and have been conveniently normalized to a simple 0(oP s ) term. As 
indicated, the leading order contribution to T% vanishes. In the presence of target mass, the 
physical structure functions in Eq. ([3]) - relating to a tensor projection on the hadronic 
W^ v - are a mixture^ of Ti. Explicitly in our case of CC charm production we have: 



F{ = T{ (14) 
f 2 C = 2 ^ (15) 



Ft = 2 S (16) 



^ = --K^^ (18) 



v 



where 



" = 1 + l^J • (19) 

Fj C = F£(x,Q 2 ) and .7-? = jF 4 c (a;, Q 2 ). These results are in agreement with Table V in 
Ref . [ 20| . They rely only on the assignment of parton light cone momentum p + related to 



the nucleon light cone momentum p + = £P/v, with p± = and £ = 77 at leading order. 

The NLO corrections 5 were first obtained in Ref. [18| and iiff'f 3 have been red- 



erived in [21]. We follow these references in notation and present the full set [i = 1, ...,5) 



5 See Ref. [ 20 1 for details on the mixing of structure functions within the tensor basis. 
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including our independently rederived fff'f . The fermionic NLO coefficient functions Hf in 
Eq. (|ToD for charm production are calculated from the subprocess W + s — *> gc and from the 
one-loop correction to W + s — > c. They are given by[] 



where P^(£) 



^1=1,2,3,5 ^2 ' A 

4 A + e 2 \ 



n ( ?(0 ln-^ + fc?&A) 
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(20) 



(21) 
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" = M 4+ 2A + Y + 2A 
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21n(l-Q-ln(l-AQ 
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(22) 



and 



K A 



[1 - A) ln(l - A) . 



(23) 



The coefficients in ( pT|) for i — 1, 2, 3, 5 are given in Table |. A misprint in Ref. ]TB| concerning 



A 2 was already corrected in [^TJ. Here, we correct a similar misprint concerning A^. 

The gluonic NLO coefficient functions Hf in Eq. ( |T0"| ) for charm production, as calculated 
from the subprocess W + g — > cs, are given by 



where pW(£) 



e 2 + (i-0' 



A) 



(24) 



In 



;i-ax 



and 



/*?(£, A) = C + C M £(1 - + C a ,i + (1 - A) £ (Cs, + A £ C 4 ,) 



(25) 



with 



Co = [21n(l-£)-ln(l-A£)-ln£] 



(26) 



6 For the ease of notation, from here on £ will be a generic variable and no longer the light-cone 
momentum of Eq. (^). To strictly match Eq. ( |i~0| ) we should have £ — > £' instead. 
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The coefficients Ck,i are given in Table [TT[ As in Ref. [^TJ, they differ from those in Ref. [18 
by counting - within dimensional regular ization - the gluonic helicity states in D [22,23 
rather than in 4 dimensions. 

The structure function F 4 is insensitive to any collinear physics up to 0(a]) and the 
corresponding coefficients H^ 4 are, therefore, scheme and scale-independent functions. They 
are given by 



4 A(l-flg[l + (l-2A)fl 
3 (1 - AO 2 



(27) 



tff =4 (£,A) = 2A£ [l-£-(l-A)£ L> 



(28) 



In the limit A — > 1 (jn c — > 0) and after an additional minimal subtraction of the collinear 
mass singularities in H? =l2 s 5 the Hf' 9 reduce to the massless MS coefficient functions:[] 



lim {//:' ,,, (C, g, A) T (1 - SjPg* In ^ j 



41(0 
41(0 



(29) 
(30) 



m 2 =Q 2 



where 4i and C^^ for i = 1,2,3 are the massless MS coefficients in Appendix III of p3" 



Extending the results in Refs. p2|-p4| to include i — 4, 5 and within the notation of Ref. |23 
we find 



4ko = 3 e 

4ko = 4ko 
4^0 = 2^(1-0 
41(0 = 4ko • 



(31) 

(32) 
(33) 
(34) 



To calculate the light quark contributions to CC v T DIS, the massless coefficient functions on 
the RHS of Eqs. fl2"9D- (|34j) are used together with the PDFs which multiply the CKM matrix 
elements |Vij^ c | 2 in an obvious modification (A — ► 1, H — > C, ...) of our Eq. (|T0|). The 



result is added to the charm production component fllOD to obtain the entire NLO structure 
function. 

It is interesting to note that the equalities (0), fl34|) guarantee that the Albright- Jarlskog 
relation (|2|) is not violated in massless QCD at NLO. Eqs. flTjj)-(|rj) and the fact that 
H$ 9 7^ Hi,' 9 for charm production manifest, of course, a violation of Eq. (Q) in the real world 
of massive target hadrons and of heavy quarks interacting through QCD. This observation 
must actually be expected to hold at all orders as we clarify in the Appendix. 



7 The choice \i = Q is made to match with the notation in [23]. Retaining a general fi ^ Q on the 
LHS of (|29|), (|30|) results in the massless coefficient functions for arbitrary scale. Obviously, this 
amounts to adding back the ln(Q 2 / fi 2 ) "splitting function times log" counter-terms. 
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III. STRUCTURE FUNCTIONS AND CROSS SECTIONS 



A. Structure Functions 

The Albright- Jarlskog relations (Eqs. (|]) and (Q)) are valid at leading order in the mass- 
less limit. Here, we show F 4 and 2xF 5 — F 2 to demonstrate violations of these relations. 
We use the CTEQ6 parton distributions [25] which include estimates of the uncertainties in 
the distributions. Thus, in Section |111 B| we will be able to quantify the error in the eval- 
uation of <7cq(v t N) that is caused by our imperfect knowledge of the PDFs. The CTEQ6 
fits are provided for \i > fj, — 1.3 GeV but can be extrapolated to lower values of /i. When 
the factorization scale squared goes below /i 2 ut = 0.5 GeV 2 , though, we choose to freeze it 
at fi 2 = /i 2 U f An alternate set of parton distribution functions used below is the Gluck, 
Reya and Vogt GRV98 set which evolves from /i 2 , = 0.4 GeV 2 with parametrizations 



provided above /i 2 ut = 0.8 GeV 2 . In principle, CTEQ6 covers the parameter space of global 
PDF-related data within conservative errors. On the other hand, the evolution with three 
quark flavours (u,d,s) of GRV98 matches better with our evaluation of <Tcc(v T N) at low 
neutrino energies where we assume light sea quarks only. In comparison, the CTEQ6 PDFs 
are evolved with a variable flavour number which is, strictly, not fully compatible with 
our approach. When we quantify the uncertainties from the PDF degrees of freedom we 
will, therefore, employ GRV98 for an independent comparison. Thereby, we can convince 
ourselves that the slight inconsistency mentioned above leads to no noticeable bias in prac- 
tice. The CTEQ6 NLO MS set of parton distribution functions is our default choice in the 
evaluations below. If not stated otherwise, the curves correspond to fi = Q. 

In Fig. [T], we show F^x, Q 2 = 2 GeV 2 ). The solid lines in the figure include target mass 
corrections of Eq. (|i~7|), while the dashed lines are with the target mass set to zero, so that 
i] —>■ x and p = 1. The leading order (LO) curve with = shows F4 = 0, even with 
charm mass corrections included. A leading order violation would occur if the initial quark 
masses were set to non-zero values. Including the target mass corrections shows the effect of 
the mixing of T± with Ti and JF 5 in Eq. (|PTD. The NLO corrections have an effect primarily 
at small-x. 

Fig. shows the violation of the second Albright- Jarskog relation, at fixed Q 2 = 2 GeV 2 . 
Even at leading order, 2xF^ — F2 7^ 0. This is due to charm quark mass corrections in 
W + s' — > c. The magnitude at small- x reflects the impact of the s' sea distribution. Target 
mass effects incorporated by Eq. (|18D are not significant. Including the NLO corrections 
makes small changes to the curve at small-x. Both of the figures show that in evaluations 
of the total charged current cross section, the naive Albright- Jarlskog relations are good 
approximations to the NLO results. This is true at low energies, where <7cc(v T N) does not 
probe small- 2 and at high energies where F45 are suppressed, anyway. 



B. Cross Sections 



The cross sections for neutrino and antineutrino charged current interactions with 
isoscalar nucleons are shown in the first two panels of Fig. [| We use our default set of 
CTEQ6 parton distribution functions with the factorization scale equal to the renormaliza- 
tion scale fi 2 , and fi 2 = Q 2 . Below fi 2 = /i 2 ut , we set /i 2 = /z 2 ut in the parton distribution 



S 



functions and in a s , but we keep the explicit Q 2 dependence in the differential cross sec- 
tion of Eq. (!]) as well as in the counter-logs ~ lnQ 2 //i 2 ut of the NLO coefficient functions. 
For Q 2 < /i 2 ut , a noticable impact of these logs is indicative of long-distance strong inter- 
action. The evaluation can, therefore, not be trusted perturbatively, whenever it becomes 
sensitive to the technical choice fi > fi cut - Most of our results below are, however, com- 
pletely insensitive to it. In Fig. |3] muon (anti-)neutrino cross sections appear with dashed 
lines, while the solid lines show the tau (anti-)neutrino cross sections. The upper curves 
show no cuts while the lower curves have W 2 = Q 2 (x _1 + 1) + M% > (1.4 GeV) 2 and 
Q 2 > 1 GeV 2 . As we show in Fig. |J for W min = 1.4 GeV, the tau neutrino CC cross 
section is fairly insensitive to the Q 2 cut of 1 GeV 2 . At, for example, E v = 20 GeV, 
crcc(Qmin = 1 GeV 2 )/<7cc(Qmin = 0) = 0.93 for u T N CC interactions. The Q 2 cut has a 
larger impact on acc^N), where o" C c(<3 m in = 1 GeV 2 )/<7cc(Qmin = 0) = 0.85 for E v = 20 
GeV, due to the fact that for nearly massless leptons, Q 2 > cuts out a larger fraction 
of the available phase space. 

The small changes in the CC cross sections with Q 2 ^ = 1 GeV 2 lead one to expect 
that non-perturbative effects at low Q 2 are unlikely to be large when one also applies the 
Wmin cut of 1.4 GeV. At low energies, without the W m i n cut, a substantial contribution to 
the cross sections comes from Q 2 < 1 GeV 2 , however, and this is precisely where the DIS 
cross section is only a rough approximation to the true cross sections with quasi-elastic and 
resonant as well as non-resonant contributions. 

The bottom panel of Fig. |3| shows the ratio of the u T N to u^N CC cross sections (solid 
lines) and the same ratio for antineutrinos. Shown are the uncut results, but the results 
with the Wmin and Q 2 cuts agree to within 3% for E v > 20 GeV. Of note is the fact that 
even at E u = 10 3 GeV, the u T N to u^N CC cross section ratio is still 5% below unity. At 
100 GeV, the ratio is 0.76. There are two reasons for the deficit in the v T CC cross section: 
the reduced phase space and the contribution of F 5 . The reduced phase space is reflected 
in the integration limits for x and y (Eqs. (|5|) and (|7J)) [27]. This is responsible for about 
half of the suppression of the v T CC cross section relative to oq,c{v 11 N'). In Eq. (|3|), the 
F 5 term appears with a minus sign, and no factor of x. Since F 5 ~ F x ~ q(x,Q 2 ), there 
is a small-x enhancement of its contribution to the cross section at high energies. The F 5 
term accounts for the rest of the suppression of the v T cross section at high energies. The 
tau mass corrections to the prefactors of Fi, F 2 and F 3 become negligible at high energies 
because the low-x rise of q(x) is tempered by factors of x or y for these structure function. 

Since mass effects for tau neutrino interactions persist to 1 TeV, it is interesting to 
compare to the case of muon neutrino CC interactions at low energies where the muon mass 
is important. Since the factor of lepton mass comes into the equations via m 2 /E v , the energy 
for muon neutrino interactions equivalent to 1 TeV for v T is E u = 3.5 GeV. At this energy, 
cr cc( zy A'A r ) including m M = 0.106 GeV (without cuts) is 2% lower than the cross section with 
m M = 0. This smaller suppression is due to the fact that at low energies, one does not get 
significant small- x contributions to F 5 . Furthermore, for u^N scattering, this energy is in 
the range where the DIS approximation to the total cross section is not reliable. 

Charm production in neutrino interactions is a small contribution. In Fig. [|, we show 
the total charged current cross section and the separate cross section for uN — > cX as a 
function of incident neutrino energy. The dashed curves are for incident u^, the solid curves 
for incident u T . At E u = 100 GeV, charm production contributes about 7% of the cross 
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section for both z/ M and v T . 

The K-factor, a comparison of the NLO to LO| charged current cross sections for incident 
tau neutrinos is shown in Fig. |^. As to be expected, NLO corrections are most significant 
near threshold where virtualities Q 2 are lower and a s is larger than at higher energies. At 
E v = 10 GeV, K =NLO/LO=1.12, reducing to K = 1 at E ~ 50 GeV for the evaluation 
with no cuts. The three curves show that with W m i n > 1.4 GeV and Q 2 > 2 GeV 2 , the 
K-factor is the same as for the uncut cross section except for E v < 10 GeV where the cuts 
improve the perturbative convergence. The K-factor with only the W m [ n cut is even slightly 
lower, confirming the expectation based on Fig. |] that perturbation theory is well behaved 
and its convergence does not have to be improved by any further Q m j n cut on top of W min . 
Overall, the K factor is reasonably close to one to indicate trustworthy perturbative NLO 
predictions, even more so when the W m in cut is imposed. 

C. Uncertainties 

In this Section we investigate a few factors that cause theoretical errors in the evaluation 
of dec- A full assessment of uncertainties will have to be based on the combination of the 
quasi-elastic and resonant channels with our DIS results. At higher energies, though, where 
DIS becomes dominant, we already can provide a good guide towards error estimates. 

An uncertainty in our evaluation of the cross section is due to the factorization scale 
dependence of the cross section at fixed order in perturbation theory. To evaluate the un- 
certainty due to the scale dependence, we have varied the factorization and renormalization 
scale (/i = fip = fiji) over a very wide range of /i 2 = 0.1 — 10 Q 2 for energies between 5 GeV 
and 100 GeV. A selection of energies is shown in Fig. |7|. We show the ratio of aQc{v r N) as 
a function of [i 2 /Q 2 to the cross section at /i 2 = Q 2 . The flat /i dependence is a reassuring 
feature of the NLO calculation as opposed to the monotonic decrease with n observed in LO, 
see Fig. ||. While we discuss here the full range of fj, 2 as plotted in Fig. [7], the perturbative 
stability observed in Fig. || suggests that one would very likely overestimate the uncertainty 
from higher orders by such a wide scale variation. This seems even more so if we look at 
the scale dependent K-factor in Fig. [| which prefers scales close to the canonical DIS choice 
fi = Q where K ~ 1. Still, scale choices are arbitrary to a large extent and we prefer to 
discuss the full picture instead of narrowing it down to some window for \i around fi = Q. At 
E v = 5 GeV, the variation is the largest for the cross section evaluated with no cuts, ranging 
between 1.16 and 0.66. As explained in Section |111 A] , for /i 2 < /i 2 ut , we always set \i = /i cu tQ- 
For small fi/Q, the ratio without cuts in {W, Q} is nearly constant as a function of fi 2 /Q 2 , 
showing the degree to which the uncut cross section comes from small /i values where /i cu t 
takes over and where the perturbative treatment is unreliable. With {W, Q} cuts applied, 
the variation over the range of scales is slightly less. The scale uncertainty at low energies 



8 The LO result is obtained by neglecting all O (a*) terms in Section || and employing the CTEQ6L 
PDFs. 

9 I.e. strictly the plot should be labeled dec (max{/j,/i cut }) /o"cc (max{Q, /J, cu t})- 
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underlies a larger uncertainty associated with the DIS approximation near threshold where 
the quasi-elastic and resonant contributions are significant. 

At higher energies, the variation in the ratio is between ~ 0.85 — 1. For E v = 20 GeV, 
the ratio in Fig. [7| is between 0.95 and 1 for fi 2 /Q 2 ~ 0.2 — 4. 

The uncertainty in the cross section due to uncertainties in the parametrization of the 
parton distribution functions (PDF) is harder to quantify than the scale dependence. One 
approach is to use different parton distribution functions than the default CTEQ6 set. As 
a comparison, the GRV98 PDF set, at E v = 20 GeV with fi = Q, yields a cross section only 
about one percent smaller than the CTEQ6 set. At E v = 5 and 10 GeV, the GRV98 cross 
sections are 10% and 4% lower than the CTEQ6 cross sections when no cuts are applied. 
The larger deviation at lower energies is due to the different high-x distributions in the two 
PDF sets. Including W > 1.4 GeV results in a deviation of 3.6% at E v = 5 GeV, even less 
for higher energies. Table |TXT| represents results for a few selected energies with a cut in W 
applied. As commented above in Section 111 A , the PDFs also differ in their treatment of 
the number of active flavors and in the value of the strong coupling constant. 

The CTEQ collaboration has also provided distributions in addition to their best fit set 
|25[. The 20 dimensional parameter space to which the PDFs are sensitive yields 40 PDF 
sets with plus/minus variations on the eigenvector directions in that space. The resulting 
error estimate on acc{^ T N) from evaluating the 40 sets is 3% at E v = 20 GeV, see Table 



I II for other values of E u . Overall, the GRV98 results lie within the uncertainty estimate 



suggested by the CTEQ6 eigenvector PDFs. We can, therefore, be confident of the absence 
of a systematic effect from the number of flavours and that the statistical uncertainties as 
encoded in the CTEQ6 sets provide a realistic PDF error estimate for acc{v T N). 

We have incorporated kinematic corrections due to including the target hadron mass 
Mn by employing the parton light cone fraction £, which equals the Nachtmann variable r\ 
|T9| at leading order in the massless quark limit. One finds that rj is much different than 
Bjorken-x at large x. For example, for Q 2 = 2 (10) GeV 2 , rj = 0.45 (0.49) at x — 0.5 and 
rj = 0.75 (0.92) at x — 1. The use of rj rather than x in the structure functions has the 
largest impact at high x and low Q 2 . 

Target mass corrections are also included via Eqs. (|I4|)-(|T8|), in which the Ti are mixed 
with target mass dependent prefactors for a given These formulae are based on the 
assignment of p+ = £P + to the light-cone momentum of the massless incident parton (p 2 = 
0) given P 2 = Mjy. The parton and nucleon are assumed to have collinear momenta, 
p± = 0. The formalism is discussed in detail in Ref. |20|| , however, other choices for the 
model of including target mass effects are possible, for example, including parton transverse 
momentum [B8[] p± ^ 0. The model of Ref. reproduces the kinematic corrections for 
the leading twist operator product expansion result [|9[ . As a consequence, our results here 
serve only as a guide to the magnitude of target mass corrections to the charged current 
cross sections. For the transverse structure function Ft of Ref. , which equals T% in the 
collinear parton approximation, the difference between the p± ^ approximation giving Ft 
and our order p° ± approximation is less than 10% over all x and Q 2 as low as 1 GeV 2 , rapidly 
falling below 5% at Q 2 =2 GeV 2 . Low energy 0{p 2 L ) effects and the target mass treatment 



in general border onto dynamical higher twists at large-x [B0-B2J. Those concepts will be 



investigated in more detail when future work will combine the DIS cross-sections with the 
non-DIS channels. 
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The effects of the target mass corrections as implemented here on F4 and 2xF§ — F%, as 
shown in Figs. H and |2|, are small. In terms of the total cross section, the lowest energies are 
most affected. In the absence of kinematic cuts, acc( h, TN)(E u = 10 GeV) for M N = is 8% 
larger than the cross section including the target mass via 77 and Eqs. (0)"(H3)- For E > 20 
GeV, the deviation is less than 5%, down to the 2% level at 50 GeV incident neutrino energy. 
When one includes W > Wmin — 1.4 GeV, the effect of target mass corrections on the v T N 
CC cross section is less than 2% for E v > 8 GeV. The reduced effect is due to the fact that the 
W min value reduces the region of integration for large x since 2M N E v y{ \ — x)> W^ in — Mjy. 



IV. CONCLUSIONS 

The NLO corrections for the v T N charged current cross section have a relatively large 
impact at low neutrino energies near threshold and less of an impact at high energies. 
We have already shown the K-factor in Fig. In Fig. [IT], we show aQc{v T N) / E u versus 
neutrino energy at NLO including target mass corrections as implemented via Eqs. (fj4[)-(|T8|) 
and charm mass correction. This is compared to the naive evaluation, neglecting masses and 
NLO corrections, where the Albright- Jarlskog relations are correct. At high energies, the 
target mass corrections are negligible. The main effect is due to the charm mass threshold. 
At lower energies, the larger QCD K-factor is compensated by the reduction in the cross 
section due to target mass and charm mass effects. 

In the theoretical evaluation of the v T N charged current cross section, the scale depen- 
dence in the PDFs and a s is a large uncertainty at low energies. The parametrization of the 
PDFs does not seem to be a large uncertainty in the evaluation of the total cross section, 
especially when one applies DIS cuts on W 2 and Q 2 . 

Target mass corrections and the importance of Q 2 < 1 GeV 2 are small at high energies, 
but they are significant in the 10-20 GeV energy range and lower, especially for the cross sec- 
tion without kinematic cuts. Our implementation of target mass corrections neglects parton 
transverse momentum (pj_). This approach makes the inclusion of NLO QCD corrections 
straightforward, however, it neglects corrections of order M^/Q 2 induced by nonzero p\ 
||28|| . As a consequence, one should view our target mass corrections as approximate. 

In any case, the DIS cross section in the energy region below E v ~ 10 — 20 GeV is 
difficult to interpret. At best, the uncut cross section is a crude approximation to the total 
cross section including resonant hadron production, e.g., A production. The cross section 
with cuts is likely a better representation of the non-resonant neutrino-nucleon interactions, 
however, the issue of avoiding double-counting in combining resonant and non-resonant 
interactions is not solved theoretically. Phenomenological approaches are being explored 
13fl for applications to v^N interactions in the few GeV region that may also be applied 



to u T N interactions at higher energies. The universality of M^/Q 2 corrections, carried over 
from copious electromagnetic interaction data to the context of weak interactions, is not 
completely clear. 

The vN cross section is an important ingredient in current and future atmospheric and 
neutrino factory experiments. Our evaluation of the NLO corrections for z/ T N CC interactions 
including charm mass corrections and an estimate of target mass effects is part of a larger 
theoretical program to understand the inelastic vN cross section over the full energy range 
relevant to experiments. 
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APPENDIX A: THE ALBRIGHT- J ARLSKOG RELATIONS 



In this Appendix we pinpoint the approximations for which the Albright- Jarlskog rela- 
tions (AJRs) (0) and (Q) hold beyond the naive parton model (where they are exact). Below, 
we will rewrite the AJRs in terms of helicity amplitudes. First, however, we can immediately 
tell from Eqs. (|3)-(0D that the mixing of {J^, T% % for M N ^ will violate both (@) and 
(§). We will, therefore, restrict the following discussion to = 0. The charm mass will 
be retained to trace its (somewhat less obvious) impact on the AJRs. In terms of helicity 
projections 



W 


= e£e%Wr 


w s 


= (4< + eg 


w q 




w± 





(IV 



(Al) 
(A2) 
(A3) 
(A4) 



with polarization vectors in terms of virtual W momentum q and longitudinal reference 
vector k 



the tensor basis Tt 



e± = -j= (0, =Fl, — i, 0) (transverse) 
n __ g M (scalar) 

can be written as: 



(A5) 
(A6) 

(A7) 



^=-{W + + W.) 

?2= + W. + 2W ) 

H= t\{W + -W.) 



(W+ + W„+ 2W - W s 



(A8) 
(A9) 
(A10) 
(All) 
(A12) 



We then see that 



r 2 



2xF' 



x 



2xTl 



xW* 



(A13) 



singles out the interference term W s between scalar and longitudinal polarization. The 
latter involves a contraction with oc g M . Now, as long as we make the single W boson 
exchange approximation, the DIS process is equivalent to an effectively abelian (electroweak) 
interaction if all quarks are massless (or of the same mass). Then, F 2 —2xF 5 = is guaranteed 
by naive gauge invariance under — > — q^ j \f—q 2 = 0^. This is a stronger statement 
than helicity conservation for massless spin- 1/2 quarks because gauge invariance holds to 
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any order in a s while helicity conservation breaks down when non-collinear NLO radiation 
generates angular momentum. A transition W + s — *> c with = m s ^ m c , however, is 
necessarily non-diagonal in flavour space and the interaction cannot be abelianized. Thus, 
naive gauge invariance does not hold and the second AJR (||]) is not protected from charm 
mass corrections. It will, however, hold at any order in a s for massless quarks. 

Worth mentioning is also that F£ = at O even for the charm production process, 
indicating that Wq + W q — W s = at LO. This should be compared to the longitudinal 
structure function F 2 — 2xF\ oc Wq which does not vanish by helicity conservation in LO 
when a massive charm quark is produced in the final state. Now, F 4 is obtained from 

W + W g -W s = -qtW, v . (A14) 

At the parton level, k is the incoming parton momentum p. F4 does, therefore, not receive 
corrections from the final state charm mass in LO as long as initial state down and strange 
masses are zero (p 2 = 0). 
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TABLES 



i 


Ai 


Bi t i 




-83,2 


1 

2 
3 
5 




K A 


X-1-Ka 
X 


1 - 4^ + «e 2 

2 - 2£ 2 - | 

-i-c 2 


3-2£-£ 2 


1 

2 
1 
2 
1 

2 



TABLE I. Coefficients for the expansion of hf in ( |2l| ) 



Cl,i C2J C3; C, 



4./ 



1 


4-4(1 - A) 


1-A£ 1 


2 


-4 


2 


8-18(l-A) 
+12(1-A)2 


1-A 1 


6A 


-12A 


3 


2(1 -A) 





-2(1-6 


2 


5 


8- 10(1 - A) 


(1-A)£ -I 


4 


-10 



TABLE II. Coefficients for the expansion of hf in (|2 



[GeV] 



A a C c 



acc(GRV98) - f j C c(CTEQ6) 



5 

10 
20 
50 
100 



5.6% 
3.3% 
2.8% 
2.4% 
2.2% 



-3.6% 
-2.5% 
-1.0% 
±0.0% 
+0.5% 



TABLE III. Propagation of PDF uncertainties into the evaluation of acc{v T N) with a cut 
W > 1.4 GeV. The second column Aucc was calculated using the CTEQ6 eigenvector PDFs 
along the master formula (3) in [25|. The third column compares GRV98 [^] with the central 
CTEQ6M set. 
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FIGURES 




X 

FIG. 1. The structure function F4 as a function of x at fixed Q 2 = 2 GeV 2 , at LO and NLO, 
using the CTEQ6 parton distribution functions. Dashed lines show the case with target mass 
M N = 0. 
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X 

FIG. 2. The structure function difference 2xF§ — F2 as a function of x at fixed Q 2 = 2 GeV 2 , 
at LO and NLO, using the CTEQ6 parton distribution functions. Dashed lines show the case with 
target mass Mjy = 0. 
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FIG. 3. Cross sections for inclusive neutrino [accivN), upper panel] and anti-neutrino 
Wcc{vN), middle panel] production of charged leptons on an isoscalar target, evaluated at NLO 
and plotted versus the energy of incident /j,- and r-flavoured (anti-)neutrinos. We show the out- 
come of a naive integration over the full kinematic range of W 2 and Q 2 along with the effect of 
imposing DIS cuts on these variables. The ratios occ{v T N) / cjqq^^N) and ucc{i> T N) / acc{^n^) 
which we show in the lower panel are insensitive to these cuts (and we show the uncut results). 
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Q^/GeV 2 

FIG. 4. Ratio acc( ( 3min)/ (T cc( ( 5min = 0) versus Q^in f° r u tN CC interactions with E u = 5, 
10, 20, 50 and 100 GeV. The solid lines are with W min = 0, the dashed lines with W min = 1.4 GeV. 
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E [GeV 



V 



FIG. 5. The charged current cross section and the separate contribution from vN — > cX for 
incident (dashed line) and v T (solid line). At E u = 100 GeV, the charm production contribution 
is about 7% of the total. 
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FIG. 6. The K-factor K =NLO/LO for tau neutrinos with no cuts (solid), W min = 1.4 GeV 
(long dash) and W min = 1.4 GeV with Q 2 > 2 GeV 2 (short dash). 
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5 ' 1 1 1 1 1 1 1 L - ' 1 1 1 1 1 1 1 

10 _1 l H 2 /Q 2 10 

FIG. 7. The ratio of <tcc(aO /^cdi 1 = Q) m NLO for z/ r N interactions, as a function of n 2 /Q 2 
for several values of E^. 
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FIG. 8. The ratio of ace (m) /^"CC (a* = Q) m LO for z/ r N interactions, as a function of /J?/Q 2 
for several values of E v . 
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FIG. 9. The K-factor K =NLO/LO versus factorization scale /i for tau neutrinos with no cuts 
(solid), W min = 1.4 GeV (long dash) and W min = 1.4 GeV with Q 2 > 2 GeV 2 (short dash). 
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FIG. 10. Our calculation of <jqq(v t N) / E u versus neutrino energy compared to its naive evalu- 
ation neglecting masses and NLO corrections. 
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